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(54) Method and apparatus for excision of narrow band interference signal from CDMA signal 



(57) A narrow band interference excision circuit for 
use in broadband digital data communication systems 
such as CDMA systems. The excision circuit is com- 
prised of a matrix of polyphase filters that divide the input 
signal into a plurality of narrow subbands. Each narrow 
subband signal is examined to determine if narrowband 
interference exists in that bin. This is done preferably by 
taking the average of the absolute amplitude of the sig- 



nals in the bin. If a signal in a bin has an absolute am- 
plitude which far exceeds the average in the bin, the en- 
tire bin signal is eliminated. A bank of polyphase syn- 
thesis filters reassembles the composite signal. An 
equalization circuit with an error predictor comprised of 
an adaptive FIR filter and a correlator is coupled to adapt 
coefficients of the filter to colored noise and generate a 
colored noise cancellation signal to remove colored 
noise from the input to the slicer. 
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Description 
Field of Use 

[0001] The invention is useful in code division multi- 
plexed digital data transmission systems and other dig- 
ital data transmission systems using carriers to elimi- 
nate or substantially reduce narrow band noise. 
[0002] In a paper by Jeffrey A. Young and James Leh- 
nert, entitled "Analysis of DFT-based Frequency Exci- 
sion Algorithms for Direct-Sequence Spread- Spectrum 
Communications" published in August 1998 in IEEE 
Transactions on Communications, Vol. 46, No. 8, p. 
1076 (hereafter the Young paper), the authors describe 
a discrete Fourier transform type frequency excision al- 
gorithm to eliminate narrow band noise from direct-se- 
quence spread-spectrum modulation. The authors note 
that processing gain limits the interference rejection ca- 
pability of unaided direct-sequence spread-spectrum 
modulation. The prior art contains numerous narrow 
band interference rejection techniques called frequency 
excision algorithms to extend the interference rejection 
capability of spread spectrum systems beyond the 
processing gain limits. This is an important ability to 
overcome severe interference situations when strong 
narrow band spurious signals are received along the 
spread spectrum signal. 

[0003] This frequency excision capability is highly de- 
sirable in, for example, the new digital data delivery sys- 
tems for delivering high bandwidth telephone service, 
video on demand and high speed internet access to sub- 
scribers on cable TV systems via the hybrid fiber-coax 
cable plant. Frequency excision also allows transmitted 
power levels to be reduced. 

[0004] The prior art adaptive notch filtering tech- 
niques described in references 1-3 of the Young paper 
are noted to be cumbersome and adapt slowly to the 
frequency of the interfering signal. Adaptive A/D conver- 
sion described in reference 4 of the Young paper works 
for a single CW interference source but cannot be easily 
generalized to multiple interference sources. Transform 
domain signal processing using SAW filters described 
in references 5-11 of the Young filter promises wide 
communication band widths and rapid adaptation to 
changing interference. However, a drawback to using 
SAW filters is that they limit the linear dynamic range 
and, hence, limit the immunity to multiple interferes. A 
class of excision algorithms called adaptive digital filter- 
ing (ADF) described in references 7, and 12-16 of the 
Young paper, is very large and consists of adaptive 
transversal filtering, classical filtering with parameter es- 
timation, lattice filtering and decision-feedback filtering. 
ADF provides a wide variety of estimation algorithms 
with varying response times, but the algorithms that 
were study by Young and his co-author were limited in 
the number of interferers that can be rejected by the 
number of delays (poles and zeroes). Young notes that 
the advantage of discrete Fourier transform (DFT) 
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based frequency excision is the ability to handle multiple 
interference sources and the ability to adapt rapidly. The 
class of DFT algorithms studied by Young are described 
in references 1 7-25 of the Young paper. DFT algorithms 

5 differ from SAW filter based techniques in that digital 
technology is used which yields greater freedom in de- 
signing the interference removal algorithm and easily 
provide high dynamic range. The number of interfering 
signals that can be removed is related to the length of 

io the DFT and may easily extend into the hundreds for a 
1 024 point DFT and notch depths on the order of 60 dB 
can be achieved. 

[0005] However DFT excision algorithms are soft- 
ware based. Because of this fact, they are too slow for 
15 many applications where data rates and traffic volume 
are very high such as in digital service delivery over ca- 
ble modem based systems. 

[0006] In a 1 994 paper by Kohri, An Interference Sup- 
pressor for CW and Narrow-Band Signals Using Filter 

20 Bank on CDMA Communications, published July 4-6, 
1994 at the University of Oulu, Oulu, Finland in the pro- 
ceedings of the IEEE ISSSTA'94, Kohri proposed a nar- 
rowband interference suppressor using a bank poly- 
phase FIR decimating (dividing) filters, a limiter and a 

25 bank of combining or interpolator filters at the front end 
of a CDMA receiver. The transfer function of the filter 
bank was a series of individual, non-overlapping trans- 
fer functions, and no error predicting equalizer was 
taught to cancel colored noise. Likewise, perfect recon- 

30 structbn filters were not taught. The .limiter was taught 
as a nonlinear amplifier. The fact that nonoverlapping 
transfer functions for the filters were used means there 
will be blind spots in the interference suppressor which 
could let narrow band interference signals through. The 

35 fact that perfect reconstruction filters are not taught, 
means that the filter banks themselves can introduce 
distortions in the spread spectrum signal which can 
cause errors in the payload data. 
[0007] Therefore, a need has arisen for a narrow band 

^o excision algorithm or machine that can deeply notch 
multiple interfering signals rapidly at the data rates of 
cable modem systems. 

Summary of the Invention 

45 

[0008] A narrow band interference excision circuit for 
use in any broadband digital data communication sys- 
tems such as CDMA or TDM A systems is disclosed 
herein. "Broadband' as the term is used herein means 

50 any transmitted signal with a broad bandwidth such as 
code division multiplexed signals or TDMA signals 
where the symbol rate is high. Basically, TDMA signals 
with a symbol rate approaching or exceeding the chip 
rate of CDMA systems has as high a bandwidth or high- 

55 er for the transmitted signals as CDMA signals. The term 
"transmitted signal" or "transmitted signals" in the claims 
is intended to include both TDMA and CDMA signals as 
well as any other broad bandwidth transmitted signal 
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that could have narrow bandwidth interfering signals 
therein. 

[0009] In the preferred embodiment, the excision cir- 
cuit is comprised of an analysis matrix of brick wall, low 
side lobe polyphase fitters and a matrix of synthesis fil- 
ters. Together, these two collections of filters implement 
a set of perfect reconstruction filters. The analysis filters 
function to divide the input signal into a plurality of nar- 
row subbands and have overlapping frequency re- 
sponses so as to eliminate blind spots in analyzing the 
entire broadband spectrum. Each narrow subband sig- 
nal is examined continuously or iterative ly to determine 
if narrowband interference exists in that bin at the time 
of each iteration. This is done preferably by taking the 
average of the absolute amplitude of the signals in the 
bin. If a signal in a bin has an absolute amplitude which 
far exceeds the average in the bin, the entire bin signal 
is eliminated until the interference signal goes away. Al- 
ternative embodiments take the average power in every 
bin or do an FFT of every bin to look for noise peaks or 
compute the variance of the signal amplitude or power 
at every frequency from the mean and, if any peak ex- 
ceeds some threshold delta value (which may be pro- 
grammable or adaptive in some embodiments), the bin 
is erased or a notch filter is programmed to take out the 
peak. Since this process is carried out iteratively on eve- 
ry bin, interfering signals which have their bandwidth in- 
crease and decrease cause as many bins as they infect 
at any particular iteration to be erased and to continue 
to be erased on subsequent iterations until the interfer- 
ence level drops below the threshold delta value for 
those bins. In other words, the process is an ongoing 
evaluation o1 every bin, and every bin that is infected 
with an interfering signal on any particular iteration will 
be erased or suppressed. 

[001 0] A bank of polyphase synthesis filters reassem- 
bles the composite signal. Polyphase filters and the No- 
ble Identity are used to enabling lowering the complexity 
of the filter structure by using decimators to lower the 
sample rate entering each filter and up converters after 
the synthesis bank to raise the sample rate back up the 
sample rate going into the decimators. In alternative em- 
bodiments, polyphase filters need not be used and more 
complex analysis and synthesis filter banks are used 
and more complex detection and cancellation circuits 
are used so as to be able to work at the higher sample 
rate. 

[0011] An equalization circuit with an error predictor 
comprised of an adaptive FIR filter and a correlator is 
coupled to adapt coefficients of the filter to colored noise 
and generate a colored noise cancellation signal to re- 
move colored noise from the input to the slicer. 

Brief Description of the Drawings 

[0012] Figure 1 is a block diagram of a typical CDMA 

receiver in which the invention finds utility. 

[0013] Figure 2 is a block diagram of the highest level 



of functionality of the narrowband excision circuit. 
[0014] Figure 3 is diagram of the frequency response 
of a single analysis filter and the subband signals within 
showing an average absolute magnitude and a narrow- 

5 band interference peak. 

[0015] Figure 4 is a diagram of the overlapping fre- 
quency responses of the filters in the analysis filter bank. 
[001 6] Figure 5 is a block diagram of one embodiment 
for a narrowband excision circuit which does not use 

10 polyphase filters and the Noble Identity. 

[001 7] Figure 6 is a block diagram of one embodiment 
for a detection and cancellation circuit. 
[0018] Figure 7 is an illustration of a combination of a 
digital filter and decimator which can be implemented 

15 using a polyphase filter. 

[001 9] Figure 8 is an illustration of a polyphase imple- 
mentation of the filter of Figure 7 with M=2. 
[0020] Figure 9 is an illustration of a polyphase imple- 
mentation of the filter of Figure 7 with M=2 after utiliza- 

20 tion of the Noble Identity. 

[0021 ] Figure 1 0 is a polyphase implementation of the 
narrowband excision circuit of Figure 5 using the Noble 
Identity, decimators and up-conve iters. 
[0022] Figure 11 is a diagram showing how commu- 

25 tator switches can be used to implement the decimators 
of Figure 10. 

[0023] Figure 12 is a diagram showing how commu- 
tator switches can be used to implement the up convert- 
ers of Figure 10. 
30 [0024] Figure 1 3 is a diagram of a prior art DFE equal- 
izer. 

[0025] Figure 14 is a diagram of a preferred form of 
equalizer for use in broadband digital communication 
systems to eliminate colored noise from the input of a 
35 slicer by using an error predictor circuit. 

[0026] Figure 15 is a diagram of a preferred form of 
error predictor for use in the equalizer of Figure 14. 

Detailed Description of the Preferred and 
40 Alternative Embodiments 

[0027] Referring to Figure 1 , there is shown a typical 
CDMA receiver in which the invention may be usefully 
employed. The filtering system of the invention is useful 

45 in any digital communication system using carriers and 
broadband waveforms where narrow bandwidth strong 
interfering signals can occur. Many spread spectrum 
systems including not only direct sequence but also 
pulsed FM chirp systems transmit their signals using 

so wide bandwidths. Those skilled in the art are advised 
that the invention is useful to eliminate narrowband in- 
terference in any broad bandwidth communication sys- 
tem, and the claims are not intended to be limited to di- 
rect-sequence spread-spectrum CDMA systems trans- 

55 mitting digital data over cable TV media, although that 
is an environment in which the invention is useful. 
[0028] In the particular application in which the as- 
signee empbys the invention, a headend or central unit 
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(hereafter CU) is coupled to one or more remote units 
or RUs by a hybrid fiber coaxial cable network of a cable 
television system. In the receiver of Figure 1 , a filter and 
down conversion circuit 10 bandpass filters the analog 
signals on channel 12 using a passband centered on 
the upstream channel frequency if the receiver is in the 
CU. The passband is centered on the downstream 
channel frequency if the receiver is in a remote unit or 
RU (frequency division multiplexing is assumed to sep- 
arate upstream and downstream digital data spread 
spectrum analog signals on the shared channel). A 
down converter in circuit 10 then converts the filtered 
analog signal back down to a baseband signal on line 
14. An analog-to-digital converter 16 then samples the 
analog signal to generate a plurality of digital samples. 
A matched filter 1 8 then digitally filters the samples back 
to their original pulse shapes. That is, the transfer func- 
tion of a matched filter is the complex conjugate of the 
signal to which it is matched. A baseband digital 
matched filter having a transfer function matched to the 
transmitted signal acts as a synchronizer and act as 
conjugate signal generators when the signal to which 
the filter is matched appears at its input. The delay line 
based matched filter is intended to recognize a particu- 
lar code sequence to which it is matched. Each delay 
segment has a delay equal to the period of the chip clock 
so that each segment contains energy corresponding to 
one chip in the sequence at any particular time. Matched 
filters and the other components of the receiver are de- 
scribed in Dixon, Spread Spectrum Systems with Com- 
mercial Applications (3rd Ed. 1994) Wiley & Sons, 
ISBNO-471 -59342-7, and Haykin, Communication Sys- 
tems (3rd Ed. 1994) Wiley & Sons, ISBN 
0-471-57178-8, both of which are hereby incorporated 
by reference. Also incorporated by reference herein are: 
Horowitz and Hill, The Art of Electronics, (Second Edi- 
tion 1984) Cambridge University Press, ISBN 
0-521-37095-7; Data and Computer Communications 
by Dr. William Stallings, Macmillan Publishing Co., New 
York (4th Ed. 1994) ISBNO-02-415441-5; Lee and 
Messerschmit, Digital Communication, (2d Ed., 1994) 
Kluwer Academic Publishers, Boston, ISBN 0 7923 
9391 0; Elliott, Handbook of Digital Signal Processing: 
Engineering Applications, (Academic Press, Inc. San 
Diego, 1987), ISBN 0-12-237075-9; Oppenheim & 
Schafer, Digital Signal Processing (1 975) Prentice Hall, 
Englewood Cliffs, N.J. ISBN 0 13 214635 5. 
[0029] After resynchronization in the matched filter a 
despreader 20 demultiplexes the spread spectrum sig- 
nal using an inverse code matrix to the code matrix that 
was used in the transmitter to spread the symbol data. 
The resulting signal on line 22 is input to an equalizer 
24 which functions to correct for the non ideal channel 
response. The corrected output signal from the equal- 
izer is error corrected by an error correction circuit using 
the redundant ECC bits added to the payload data in the 
transmitter. A framer circuit 28 re-assembles the data 
into frames or packets of the type used by the upper 
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layer protocols such as by constructing MPEG2 packets 
from the error corrected data on line 30. The resulting 
reframed data is output to the upper layer protocols for 
further processing on line 32. 

5 [0030] To employ a narrowband excision circuit in the 
receiver of Figure 1 , the invention employs a filter circuit 
like that shown in Figure 2 on bus 21 coupling the 
matched filter 1 8 to the despreader 20. The narrow band 
interference excision filter of Figure 2 is comprised of a 

10 bank of narrow bandpass filters 34 each of which has a 
different center frequency and a passband which over- 
laps slightly with the passbands of its neighboring filters 
on either side so as to prevent blind spots. 
[0031] The function of the bank of bandpass filters 34 

15 js to provide a plurality of output signals each of which 
represents the energy of the incoming signal in a narrow 
frequency range corresponding to the passband limits 
of that particular filter. This allows narrow bandwidth in- 
terference signals to be isolated in individual frequency 

20 bins for analysis and excision. In the preferred embod- 
iment, 256 individual subband filters with overlapping 
frequency responses are used to cover a 4 mHz wide 
CDMA signal bandwidth. The number of filters used is 
a tradeoff depending upon the bandwidth to be covered, 

25 the amount of hardware complexity and associated cost 
that can be tolerated and the typical bandwidth of the 
interfering signals most commonly encountered. Ideally, 
the frequency response selected for each filter will be 
such that the subbands are not much wider than the 

30 most frequently encountered interference signal so as 
to prevent loss of useful signal information when an en- 
tire bin is suppressed because of the presence of an 
interference signal. However, use of fewer analysis fil- 
ters with wider bandwidths is acceptable because of a 

35 loose coupling between the bandwidth of the individual 
subbands and loss of payload data when an entire sub- 
band is suppressed. This result follows because the en- 
ergy of the payload data is spread out among all bins so 
the loss of any one or a few of them is not terribly dam- 

40 aging to recovery of data in the receiver. 

[0032] After separating the input waveform into sub- 
band energy components in a plurality of narrow fre- 
quency ranges (hereafter referred to as bins), a detec- 
tion and cancellation circuit examines each bin to deter- 

45 mine if a narrow bandwidth interference signal is present 
in that bin. If so, the amplitude of all the signals in the 
bin is reduced to zero in the preferred embodiment. In 
other embodiments, the amplitude of all frequency com- 
ponents in the bin including is attenuated to some small 

50 amplitude. This reduces the amplitude of the narrow 
band interference source to a level which is low enough 
to insure that it will not cause enough interference to 
exceed the error detection and correction capability of 
the ECC bits added to the payload data 

55 [0033] In the preferred embodiment, the method car- 
ried out by the detection and cancellation circuit of de- 
termining if a narrow band interference signal is present 
in a bin is by taking an average of the absolute ampli- 
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tudes of the waveform at each frequency in the bin. The 
absolute ampiitude of the waveform at each frequency 
in the bin is then compared to this average. If at some 
frequency, a peak is found with an amplitude which ex- 
ceeds the average by more than some predetermined 
amount (a programmable value in some embodiments), 
then the bin is deemed to have an interference signal 
present. This notion works in a CDMA system because 
spread spectrum systems are very wide bandwidth sys- 
tems where the energy of the payload data is spread out 
over the wide bandwidth. This tends to cause the trans- 
mitted signal to have an average amplitude quite similar 
to white noise. This means that every bin will have a 
signal in it that does not vary wildly in amplitude, and 
the average amplitude of the signals in all bins will be 
almost the same. The higher the number of spreading 
codes that are used, the more true this becomes. Thus, 
a narrow bandwidth interference signal will stand out like 
the proverbial "sore thumb". 

[0034] The calculation of the average absolute ampli- 
tude of the signal at each of the frequencies in each bin 
is an adaptive process which is ongoing in realtime. The 
same is true for each of the alternative methods de- 
scribed herein of determining when there is an interfer- 
ence signal present 

[0035] if an interference signal is found in a bin, an 
attenuator or other suitable circuitry in detection and 
cancellation circuit 36 causes the amplitude of the wave- 
form at all frequencies in this bin to be attenuated to zero 
or some other small value. Figure 3 illustrates this con- 
cept. In alternative embodiments, instead of attenuating 
every frequency in the bin, a narrow bandwidth band- 
stop or notch filter may be used. If a digitally tunable 
notch filter is associated with each detection and can- 
cellation circuit assigned to a particular bin, the compar- 
ison process, upon locating an interference peak in the 
bin, can send a digital parameter to the notch filter to 
center its notch on the center frequency of the interfer- 
ence signal. In some embodiments, the comparison 
process determines how wide the bandwidth of the in- 
terference signal is as well as its center frequency, and 
sends digital parameters to the notch filter to move its 
center frequency to the center frequency of the interfer- 
ence and which also alter the frequency response of the 
notch filter to adjust the bandwidth of the stop band or 
notch to correspond with the bandwidth of the interfer- 
ence signal. 

[0036] The notch filters in detection and cancellation 
circuit 36 may be in an array having one filter which is 
associated with each bin. Since it is unlikely that every 
bin will have an interference signal therein, it is also pos- 
sible to save on hardware complexity and use an array 
of notch filters which is not as large as the number of 
bins. Individual filters in this array may be selectively 
connected to the individual conductors of bus 38 to ex- 
cise the interference signals out of the signals on the 
individual lines of bus 3B in response to control data from 
the circuitry that does the detection process. 



[0037] In an alternative embodiment, an array of 
notch filters may be incorporated in the synthesis bank 
of filters so as to be selectively coupled to filter the com- 
bined output signal on line 42 in response to control data 

s from the detection circuitry in block 36. The switching 
circuitry must be such that if two or more interfering sig- 
nals are detected by the detection circuit 36, a like 
number of notch filters may be tuned to have their notch- 
es centered on the center frequency of the interference 

10 signals and then connected in series with bus 42 to filter 
each interference signal out of the combined signal. 
[0038] Referring to Figure 3, there is shown a frequen- 
cy versus amplitude plot of the signal in one bin which 
happens to have a narrow bandwidth interference signal 

15 present. Line 35 represents the shape of the filter pass- 
band that creates the bin. This skirt shape is defined by 
the transfer function of the filter. The waveform of the 
spread spectrum signal is represented by signal 37. The 
signal 41 represents a narrowband interference signal. 

20 The dashed line 39 represents the average absolute 
amplitude of signals 41 and 37 combined by superposi- 
tion. To detect the presence of the interference signal 
41 , the detection and cancellation circuit compares the 
amplitude of the combined signals 41 and 37 at every 

25 individual frequency in the bin and if the absolute ampli- 
tude of the signal waveform at any frequency in the bin 
exceeds a predetermined amount, the bin is deemed to 
contain an interfering signal. 

[0039] The process carried out by the detection and 

30 cancellation circuit is carried out in every bin or subband 
created by the bank of filters 34, so multiple interfering 
signals can be removed simultaneously 
[0040] In alternative embodiments, other ways of de- 
tecting the presence of an interfering signal in a bin may 

35 also be used. For example, the average power of the 
waveform may be calculated by squaring the amplitude 
of the amplitude at every frequency, summing the 
squares and dividing by the number of frequencies. The 
power of the waveform at every frequency is then com- 

40 pared to the average power for the bin, and if any fre- 
quency has a peak having a power level which is greater 
than a predetermined amount above the average, that 
bin is deemed to have an interfering signal. 
[0041] Another way of detecting an interfering signal 

45 in a bin is to do a Fourier transform on the samples that 
define the waveform in every bin and then look at the 
resulting frequency components for noise peaks which 
exceed the amplitude of the other frequency compo- 
nents in the bin by a predetermined amount. 

50 [0042] Another way of detecting the presence of an 
interfering signal is to calculate the average absolute 
amplitude in the bin and then compute the variance of 
the amplitude of the combined signal at every frequency 
from the average. If any amplitude exceeds a predeter- 

55 mined amount of acceptable variance, the bin is 
deemed to have an interference signal present. 
[0043] Another way of detecting the presence of an 
interfering signal is to compute the average amplitude 
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or power in each bin and compare the average to the 
averages of each of the other bins. Any bin with an in- 
terfering signal will have an average which is substan- 
tially higher than the averages of bins which do not have 
interfering signals present. If any bin has an average 
that is higher than the other averages by a predeter- 
mined (possibly programmable) amount, that bin is 
deemed to have an interfering signal present. 
[0044] After the interfering signal or signals have been 
removed, the resulting signals in the subbands on the 
buses are input to a synthesis bank of filters 40. The 
function of the synthesis bank of filters is to put all the 
component signals in the time domain on buses 38 back 
together as a single composite signal on bus 42. 
[0045] The combined system comprised of filter bank 
34, detection and cancellation circuit 36 and synthesis 
bank of filters is called a quadrature mirror filter (QMF) 
bank comprised of known digital FIR filters for the anal- 
ysis bank 34 and the synthesis bank 40. Not just any 
FIR filters will work however for these two banks of fil- 
ters. It is important that the coefficients of these filters 
be selected such that the combination of analysis and 
synthesis filters in the filter banks fall within the subclass 
of FIR filters called perfect reconstruction filters. This 
means that if the analysis bank and the synthesis bank 
are connected back to back (with the detection and can- 
cellation circuit eliminated) the output signal generated 
from the samples output by the synthesis bank 40 on 
bus 42 would be identical to the input signal defined by 
the samples on bus 31 with some constant delay and 
constant gain present. Perfect reconstruction filters are 
used for best performance in the preferred embodiment 
of the invention. Although other filters could be used 1or 
the analysis and synthesis banks, the performance 
would not be optimal and may even be unacceptable. 
[0046] However, the "perfect reconstruction filters" 
that are preferred are not actually perfect. Because a 
finite resolution (the number of filters M in the analysis 
bank 34 is not infinite) and because the filter response 
of the analysis and synthesis filters is not perfect, even 
after excision and reconstruction, there will still be some 
residual interference on bus 42. This interference can 
be modelled as colored noise. This colored noise will be 
eliminated by an error prediction circuit that will be de- 
scribed below. Since the main energy of the narrow 
band interference signal was eliminated by the QMF fil- 
ter system, the assumption can be made that the error 
prediction system works with no decision errors. 
[0047] Figure 4 is a graph showing a typical set of fre- 
quency responses for the analysis filters of filter bank 
34. Although the preferred embodiment uses overlap- 
ping frequency responses defining the subbands to 
avoid blind spots, other embodiments can use margin- 
ally overlapping or non-overlapping frequency respons- 
es or frequency responses that overlap more than is 
shown. The number of analysis filters in the analysis 
bank 34 determines the bandwidth of the frequency re- 
sponse of each filter. 



[0048] Figure 5 shows a more detailed block diagram 
of the QMF filter structure of Figure 2. The analysis bank 
of filters 34 is comprised of a plurality of M narrow pass- 
band "brickwall" bandpass filters with very low side 

5 lobes designated hl[n] where I increments from 0 to M- 
1. Each of these filters (when combined with its paired 
synthesis filter less the intermediary detection and can- 
cellation circuit) is a filter from the known subgenus of 
FIR filters known as "perfect reconstruction filters" or 

io "near-perfect reconstruction filters" and is selected to 
have low side lobes in its frequency response. Many 
■brickwall filters" with sharp skirt rolloff and narrow pass- 
bands have high side lobes. Use of filters with high side 
lobes is not desirable for an excision circuit since when 

15 a bin is found to have a narrowband interference source 
therein, it is desired to remove all energy from the inter- 
ference source. When brickwall filters are used with high 
side lobes, energy from the interference source is not 
limited to the bin in which the interference source was 

20 found but also spills out into adjacent bins. Thus, an ex- 
cision circuit would have to not only erase the bin in 
which the energy was found but also erase all adjacent 
bins which overlap with a sidelobe of the filter whose 
main passband passed the interference signal. This is 

25 too high a penalty in lost payload data. Thus, brickwall 
filters with high side lobes are the preferred form of filter 
for the analysis filters. 

[0049] Any "perfect reconstruction filter" or "near per- 
fect reconstruction" filter from the prior art which is de- 

30 signed with coefficients to define a narrow "brickwall" 
passband frequency response with low side lobes will 
suffice for the analysis and synthesis filters. Low side 
lobes, means that the frequency response has a pass- 
band like that shown in the idealized frequency re- 

35 sponse curve 35 in Figure 3 with little or no side lobe 
activity to let signals or components having frequencies 
below F1 or above F2 through to the output. Preferably, 
each of these filters 44, 46, 48 etc. is from the genus of 
FIR or digital finite impulse response filters, although 

40 digital MR filters could also be used, and, if the receiver 
front end was analog, analog SAW filters could also be 
used. 

[0050] Decimators 50, 52 and 54 etc. function to lower 
the sample rate. The sample rate of the wide band signal 

45 on bus 31 has to be high enough to satisfy the Nyquist 
criterion to prevent aliasing. However, the FIR analysis 
filters 44 etc. modify the samples in accordance with 
their transfer function to output samples which define a 
narrow frequency band. The sample rate is maintained 

50 however, and it does not have to be that high when work- 
ing at with the narrow bandwidth subband signals. 
Therefore, each decimator lowers the sample rate by 
eliminating, for example, every other sample. This al- 
lows the circuitry to be simpler and less expensive to 

55 build. The decimators decimate by a factor of M thereby 
taking only every Mth sample and ignoring the rest. 
There is no synchronization between the decimators as 
to which Mth samples in the sequences they receive 
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they each take, i.e., the Mth samples taken by the se- 
quence of decimators do not form a sequence of se- 
quential samples. 

[0051] Figure 6 is a diagram of one embodiment for 
the detection and cancellation circuits 56, 58 and 60 that 
comprise block 36 in Figure 2. The narrowband samples 
arrive on bus 61 and are coupled to a detection circuit 
64 that computes the average absolute amplitude or 
other criteria used or does an FFT on the samples from 
the bin. The detection circuit controls the position of a 
switch 66 via a control signal symbolized by dashed ar- 
row 68. If the detection circuit detects no interference 
signal present, it leaves switch 66 in the position shown 
so the samples pass unattenuated through the circuit 
and are output on bus 63. If the detection circuit deter- 
mines that an interference signal is present, signal 68 is 
altered to cause switch 66 to connect bus 63 to ground- 
ed terminal 70 thereby completely eliminating the sam- 
ples on bus 61 from the group of samples that will be 
converted by synthesis bank 40 back into a single signal 
represented by a single group of samples. An alternative 
embodiment for the cancellation portion of the detection 
and cancellation circuits would be a switching circuit 
which imposes an adaptable notch fitter in the signal 
path when a narrowband interference source is found 
and sets the coefficients of the notch filter to substan- 
tially match the center frequency of the notch filter with 
the center frequency of the narrowband interference sig- 
nal. In some embodiments, the notch filter could have a 
fixed bandwidth and attenuation value and in other em- 
bodiments, the bandwidth and/or attenuation value 
could be programmable or adaptable. 
[0052] The combination ol FIR filters and decimators 
to the left of the detection and cancellation circuits in 
Figure 5 can be implemented as known polyphase filters 
which have lower complexity. This is made possible by 
lowering the sample rate using the down-converters 50, 
52 and 54 etc. and raising it back up using up-converters 
72, 74 and 76 etc. Further, the lower sample rate also 
reduces the complexity of the detection and correction 
circuits 56, 58 and 60 etc. since each can work at a sam- 
ple rate lower by a factor of M than the sample rate on 
bus 31. For example, the detection circuit 64 will have 
to far fewer mathematical operations and comparisons 
to calculate the average amplitude or power or do an 
FFT where there are fewer samples per second. This 
makes the algorithm simpler and reduces the perform- 
ance requirements of the processor doing the calcula- 
tions to get all the mathematical manipulations that need 
to be done in real time done within the time available to 
do them. It is simply easier to do all this mathematical 
manipulation at the lower sample rate especially since 
no information is being lost by lowering the sample rate. 
[0053] In alternative embodiments, where the higher 
complexity and higher speed requirements on the 
processing circuitry to do detection and cancellation can 
be tolerated, the down-converters and up-converters 
can be eliminated. 



[0054] The sample outputs from the detection and 
cancellation circuits are input to up-converter circuits of 
which up-converters 72, 74 and 76 are typical. The func- 
tion of these up-converters is to raise the sample rate 
5 back up to the rate of bus 31 by inserting zeroes where 
the omitted samples were on each of buses 63, 65 and 
67 etc. 

[0055] The sample streams output by each of the up- 
converters 72, 74 and 76 etc. are input to another filter 

10 g|[n] forming part of the synthesis filter bank, of which 
filters 78, 80 and 82 etc. are typical. Each of these gl[n] 
filters is preferably an FIR filter having a transfer function 
which is the inverse of the transfer function of its paired 
filter hl[n] in the "perfect reconstruction filter". If the two 

15 transfer functions were convolved in the time domain, 
the result would be one. The reason the system of Fig- 
ure 2 is called a Quadrature Mirror Filter is because the 
g filters each have transfer functions in the time domain 
which are each a mirror of the corresponding h filters. 

20 The inverse transfer functions between the h and g fil- 
ters are what makes the pair of filters a perfect recon- 
struction filter. 

[0056] The theory of polyphase filters is that any dig- 
ital filter regardless of the number of its coefficients can 

25 be broken down into a plurality of subfilters which, when 
their outputs are combined, yield the same signal as the 
original filter would have yielded. The equations defining 
the nature of polyphase filters and the Noble identity are 
well known, but are included here for completeness. 

30 There are two functions H(z) and G,(z) in the frequency 
domain which define polyphase filters. They are related 
by the following relationship: H(z) = the summation of z- 
l*G,(z M ) for l=0 to l= M-1 where z is the frequency vari- 
able in the Irequency domain and M is the number of 

35 subband filters to be used. There is a relationship be- 
tween the frequency domain function G|(z M ) and the 
time domain function gl[n] as follows: Gj(z) = the sum- 
mation for n incrementing from negative infinity to plus 
infinity of g|[n]*z _ n with gl[n] = to h[Mn =1] .... with I rang- 

40 ing from less than or equal to less than or equal to M-1 . 
[0057] Figure 7 represents an implementation exam- 
ple for a decimating filter with M = 2. In Figure 7, block 
80 represents a standard filter having a frequency do- 
main representation of its transfer function equal to H 

45 (z). Block 82 represents a decimation by 2 where every 
other sample is eliminated. 

[0058] The polyphase representation of the decimat- 
ing filter of Figure 7 with M = 2 is shown in Figure 8. In 
Figure 8, two filters 84 and B6 with transfer f unctbns in 

50 the frequency domain of G 0 (z 2 ) and G 1 (z 2 ), respectively, 
each have their outputs decimated by a factor of 2 by 
decimators 88 and 90, and the outputs are combined on 
bus 92. Note that filter 84 works on the nth sample and 
filter 86 works on the n-1 sample with both filters working 

55 at the high sample rate of bus 94. 

[0059] By the Noble Identity, the structure of Figure 8 
can be reversed so that decimators 88 and 90 are cou- 
pled to bus 94 and serve to lower the sample rate by a 
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factor of 2 with their output sample streams coupled to 
the data inputs of two filters 96 and 98 having transfer 
functions in the frequency domain of G 0 (z) and G 1 (z), 
respectively, which are related by the equations above 
to the transfer functions H(z) of filter 80 in Figure 7. This 
allows filters 96 and 98 to be simpler since they can work 
at half the sample rate of bus 94. If filter 80 has 10 co- 
efficients, the filters 98 and 96 will each have 5 coeffi- 
cients, which taken together, are the same as the 10 co- 
efficients of filter 80 in Figure 7. Likewise, if M is equal 
to 10 and filter 80 has 1 0 coefficients, each of the 10 G { 
(z) filters in the structure like that of Figure 9 which would 
result from using the Noble Identity would each have 
one coefficient. Where 10 coefficients are used in filter 
80 and M = 2, filter 96 would have the odd numbered 
coefficients and filter 98 would have the even numbered 
coefficients. 

[0060] The concept of using polyphase filters and the 
Noble Identity can be extended to any number M, e.g., 
M = 3, 4 etc. Typically, 256 filters are used thereby great- 
ly simplifying filter construction and allowing filtering in 
real time of very high data rates to cancel narrowband 
noise. It is not necessary to use the noble identity, but it 
is preferred since it allows the digital filters to work at a 
much lower sample rate. Use of the polyphase filter 
technique in conjunction with application of the noble 
identity allows a simpler physical structure that can be 
implemented by less complex and expensive filters to 
implement the h,(n) and the g,(n) filters in Figure 5. By 
using polyphase filter implementation and the Noble 
Identity, the structure of Figure 5 can be implemented 
as the structure shown in Figure 10. 
[0061] In Figure 10, a series of decimators 100, 102 
and 1 04 etc. down convert the sample rate on bus 31 
by a factor of M. Each resulting stream is input to a sep- 
arate data input of an analysis polyphase filter matrix 
106. These separate data inputs are represented by 
buses 105 and 107 etc. This matrix is comprised of a 
collection of type 1 polyphase filters, each polyphase fil- 
ter doing the job of a single h,(n) subband filter in Figure 
5 and having a structure such as is shown in Figure 9 
but expanded to whatever value is selected for M. Each 
analysis filter H^z) can be represented by a type 1 poly- 
phase filter as the summation from I =0 to l= M-1 of z- 
IE w (z M ). This defines a matrix of filters defined by the 
matrix comprised of a plurality of rows, the first row start- 
ing with E 00 (z) and ending with E 0 M . n (z) and the first 
column starting with Eoo(z) and ending with E M-1 0 (z), 
and the last row starting with E M _-, 0 (z) and ending with 
: -m , M-i ( z )- Likewise, each synthesis filter F k (z) can be 
presented by a type 2 polyphase filter represented as 
v a summation from I = 0 to I = M-1 of z- < M " 1 - r )R kJ (z M ). 
The matrix of polyphase filters defined by this relation- 
ship to implement the bank of synthesis filters is com- 
prised of comprised of a plurality of rows, the first row 
starting with R 0 o( z ) and ending with R 0 ,m-i( z ) and tne 
first column starting with Roo(z) and ending with R M . 1>0 
(z), and the last row starting with Rm-i,o( z ) and ending 



with R M -i.m-i(z)- 

[0062] Each polyphase filter has a separate output, 
represented by buses 108, 110 etc. These outputs are 
coupled to separate inputs of a detection and cancella- 

5 tion matrix 112 which is comprised of a collection of in- 
dividual detection and cancellation circuits such as 
those shown in Figure 6. Each detection and cancella- 
tion circuit has its own output, represented by buses 
114, 116 etc. These outputs are coupled to individual 

10 data inputs of a synthesis polyphase filter matrix com- 
prised of a collection of type 2 polyphase filters, each 
type 2 polyphase filter doing the job of one of the gl(n) 
filters in Figure 5. The outputs of the synthesis poly- 
phase filters such as buses 120 and 122 are coupled to 

15 data inputs of up-converters represented by blocks 1 24, 
126 etc. The output sample streams of the up-convert- 
ers are combined on bus 1 28 for output to the rest of the 
receiver circuitry for despreading etc. 
[0063] The decimators and up-converter structures of 

20 Figure 1 0 can be represented by the commutator struc- 
tures of Figures 11 and 12, respectively Figure 11 illus- 
trates one way of implementing the decimators 1 00, 1 02 
through 104 in Figure 10 using a commutator switch rep- 
resented by movable vector 130. The input sample 

25 stream on bus 31 enters the input of the switch. The in- 
dividual data inputs 105, 107 etc. of the analysis filter 
matrix 106 are represented by the lines to the right of 
the movable vector 130. The movable vector 130 rep- 
resents a switch arm that moves selectively to couple 

30 bus 31 to each one of the data inputs 105, 107 through 
1 09 in a sequence so that sequential digital samples are 
coupled to sequential data inputs. Of course, a mechan- 
ical switch arm is not actually used because it would not 
able to move fast enough, so switch arm 1 30 is symbolic 

35 of any form of electronic switching circuit that functions 
in the way described herein. 

[0064] The commutator switch arm is operated so as 
to sequentially connect samples to sequential inputs. 
For example, when the -(M-1) sample (the sample that 

40 arrived at the switch input M-1 samples in time before 
the 0th sample) arrives at the input of switch 1 30, switch 
130 is operated to couple bus 31 to the Mth data input 
109. When the next sample arrives, the switch 130 is 
operated to couple the sample to the next data input up 

45 (not shown) in the sequence. When the -1 sample ar- 
rives, (the sample just before the 0th sample), the switch 
130 is operated to connect bus 31 to input 107. When 
the 0th sample is present at the input to switch 1 30, the 
switch is operated to connect bus 31 to input 105. Next, 

50 when the 1st sample arrives, the switch 130 recouples 
bus 31 to data input 1 09 such that the 1 st sample is input 
to the analysis switch matrix on input 109. This process 
is continued until the Mth sample arrives and is coupled 
to input 105 and so on. 

55 [0065] The commutator switch 1 32 used to implement 
the up-converters 1 24, 1 26 through 1 29 shown in Figure 
12 works in the same way as switch 130 to reassemble 
the distributed sample sequences on outputs 120, 122 
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through 123 into a single sample sequence on bus 128 
with a sample rate the same as the sample rate on bus 
31 . The switches 132 and 130 operate asynchronously 
however and there is no coordination between them. 
The notation used to define the symbols on outputs 1 20, 
122 and 123 also is different than the notation used to 
define the symbols in Figure 11 to indicate this asyn- 
chronous nature of operation. In Figure 12, the symbol 
on the left of the series is the earliest in time. In Figure 
1 1 , the symbol identified as -(M-1 ) is the first symbol re- 
ceived. In Figure 12, the symbol identified as 0 is the 
first symbol received. The switch arm 1 32 starts out in 
a position to couple output 1 23 to bus 1 28 when the first 
symbol (symbol 0) is received. The next symbol in time 
to be coupled to bus 1 28 is symbol 1 (not shown) which 
arrives on the next line up from line 123 (also not 
shown). The switch 132 continues to work its way up 
coupling each output line to bus 128 during successive 
symbol times which the M-1 symbol being the last sym- 
bol to be coupled to bus 128 before switch 1 32 re-con- 
nects line 123 to bus 128 to transmit symbol M thereto. 
[0066] There are subclasses of perfect reconstructed 
filters such as "cosine modulated" and "lapped orthog- 
onal". The specific perfect reconstruction filter that was 
selected for the preferred embodiment was a bank of 
cosine modulated filters with 256 filters. Any of the other 
subclasses of perfect reconstruction filters will also work 
as will a different number of analysis filters. The actual 
filter coefficients will be defined by the fact that 256 co- 
sine modulated filters are used, each filter having 312 
coefficients. The number of coefficients is also not crit- 
ical, but 312 coefficients is believed now to be best. Any 
MR perfect reconstruction filler will also suffice to prac- 
tice the invention. The filter frequency response charac- 
teristics should overlap however to prevent blind spots 
regardless of whether FIR or MR perfect reconstruction 
filters are used and regardless of which subclass is se- 
lected regardless of the number of filters or coefficients 
used. 

EQUALIZATION 

[0067] Because the number of analysis filters used is 
not infinite and because the analysis and reconstruction 
filters are not perfect, there will be some residual noise 
at the output of the excision circuitry on bus 1 28 in Figure 
10. If the input noise is CW, the residual output noise 
will be C W, and if the input noise is narrowband, the re- 
sidual output noise will be narrowband. This residual 
noise can be modelled as colored noise, i.e., noise 
which has some correlation. 

[0068] To understand the inventive equalizer, consid- 
er the prior art DFE equalizer structure shown in Figure 
1 3. This structure is comprised of a feed forward equal- 
izer 200 receiving input data, processing it and output- 
ting revised data driving one input of a summer 202. The 
difference input 210 of the summer 202 is driven by a 
feedback equalizer 206 which takes its input from the 



output of a slicer 204 which has an input driven by the 
output of the summer. 

[0069] The basic concept upon which the inventive 
equalizer is built is to "whiten" any colored noise on bus 

5 212 by making modifications to the prior art equalizer. 
Colored noise will result on bus 212 il narrowband inter- 
ference is present or even if narrowband interference is 
not present but FFE 200 is not perfect. 
[0070] The general scheme of modification to the prior 

io art DFE equalizer to whiten the noise on bus 212 is to 
add an error prediction circuit which functions to predict 
the next colored noise interval from the previous colored 
noise interval. The arrangement for the preferred form 
of equalizer is shown in Figure 14. Data enters the 

15 equalizer 24 on line 22 and is processed by a conven- 
tional feed forward equalizer 200. The output of the feed 
forward equalizer on line 212 is coupled to the nonin- 
verting input of a summer 214. The inverting input 210 
of summer 214 is coupled to the output of conventional 

20 feedback equalizer 206. The input 208 of the feedback 
equalizer is coupled to the output of a conventional slicer 
204. 

[0071] The output of the summer 21 4 on bus 220 has 
digital samples that define a signal which has both pay- 
25 load data plus white noise plus colored noise thereon. 
It is the function of the equalizer 24 to "whiten" the color- 
ed noise so that it does not adversely affect detection of 
each symbol by the slicer. 

[0072] The way the equalizer whitens the colored 

30 noise is as follows. First, the input samples on bus 22 
are processed by the feed forward equalizer 200 in a 
conventional manner to compensate for some impair- 
ments in the channel. The output of the feed forward 
equalizer on bus 21 2 is summed with the output on bus 

35 21 0 of a feedback equalizer 206 in summer 21 4 in a con- 
ventional manner. The output of summer 21 4 on bus 220 
is summed with correction data on bus 222 in summer 
224, and the output data is input to conventional slicer 
204. The slicer outputs the actual data transmitted plus 

40 error correction bits that were transmitted with the data 
on bus 23. The data on bus 23 is applied to the input 
208 ol the feedback equalizer 206 which processes it in 
conventional fashion to generate the feedback equali- 
zation data on bus 210. 

45 [0073] The whitening of the colored noise starts by 
coupling the signals on bus 220 to the noninverting input 
of a summer 226. The signals on bus 220 will include 
data, white noise and colored noise. Simultaneously, the 
output data from the slicer on bus 23, which is data only, 

so is coupled to the inverting input of the summer 226. The 
summer 226 functions to subtract the data only signal 
on bus 23 from the combined data, white noise and 
colored noise signals on bus 220. The resulting signals 
on bus 228 are white noise and colored noise only. 

55 These signals are input to an error predictor circuit 230. 
[0074] The error predictor circuit is a combination of 
a correlator and an adaptive digital FIR filter which has 
its coefficients adapted by the process of correlation of 
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the colored noise. Figure 1 5 is a block diagram of the 
error predictor 230. The function of the error predictor 
circuit is to adapt the coefficients of its FIR filter in ac- 
cordance with the correlation of the colored noise so as 
to predict the colored noise. The resulting data output 
on bus 222 is applied to the inverting input of summer 
224 to subtract out the colored noise from the signals 
on bus 220 so as to eliminate colored noise at the input 
232 of the sheer. In Figure 15, the adaptive FIR fitter is 
shown at 234. Its input is coupled to receive the white 
noise and colored noise signals on bus 228. The coef- 
ficients of the FIR filter are adjusted by data on bus 236 
output by a correlator 238. 

[0075] Referring again to Figure 1 4, a summer 240 
receives the data plus white noise plus some colored 
noise on bus 232 (assuming the colored noise has not 
yet been canceled). The inverting input of summer 240 
is coupled to bus 23 to receive the data. The summer 
subtracts this data from the data, white noise and color- 
ed noise on bus 232 to output a signal on bus 242 which 
is white noise and any residual colored noise. The sig- 
nals on bus 242 are input to the correlator 238 in Figure 
1 5. If there is any residual colored noise on bus 242, the 
correlator will generate a nonzero correlation signal on 
bus 236 thereby adapting the coefficients of the filter 234 
in a direction to generate an output signal on bus 222 
which causes summer 224 to subtract out at least some 
of the colored noise on bus 220. The adaption process 
continues changing the coefficients of the filter 234 until 
a signal on bus 222 causes cancellation of all the color- 
ed noise at the input 232 of the slicer. This is the desired 
state of convergence since colored noise entering the 
slicer can cause errors in determining which symbol or 
constellation point was sent during each symbol time. 
Thus, when the convergence process has been com- 
pleted, the signals on bus 232 will be data plus white 
noise only, and the signals on bus 242 will be white noise 
only. Since white noise has no convergence value be- 
tween samples, the output of the correlator on bus 236 
will be zero so no further adaptation of the filter coeffi- 
cients occurs between samples. This causes the filter 
output on bus 222 to remain at a stable sample value 
adequate to remove colored noise from the signals on 
bus 220. When the colored noise correlation properties 
of the signals on bus 220 change from what they were 
to cause the convergence, the convergence process 
starts again until the new or altered colored noise sig- 
nals are removed from the input signals to the slicer by 
summer 224. 

[0076] Although the invention has been described in 
terms of the preferred and alternative embodiments dis- 
closed herein, those skilled in the art will appreciate nu- 
merous modifications and improvements without de- 
parting from the spirit and scope of the invention. All 
such modifications and improvements are intended to 
be included within the scope of the claims appended 
hereto. 
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5. 



6. 



A method for removing narrowband interference 
signals (41) from transmitted signals (37), charac- 
terized in that It comprises the steps of: 

dividing (34) the transmitted signal up into a 
plurality of narrow subband signals with suffi- 
ciently overlapping frequency limits to avoid 
blind spots or spots of limited sensitivity in the 
broadband spectrum to be analyzed for the 
presence of narrowband noise; 
analyzing (36) the signals in each subband to 
determine if a narrowband interference source 
is present in said subband; 
removing or attenuating (36) all narrowband in- 
terference sources from all subbands in which 
they are found; and 

reconstructing (40) the transmitted signal from 
the individual subband signals after any nar- 
rowband interference signals have been re- 
moved. 

The method of claim 1 wherein the step of dividing 
(34) the transmitted signal into subbands comprises 
the step of passing the transmitted signal through a 
bank (106) of polyphase filters (84, 86, 96, 98) with 
overlapping frequency responses and different 
center frequencies. 

The method of claim 1 wherein the step of dividing 
(34) the transmitted signal into subbands comprises 
the step of passing the transmitted signal through 
the first filter bank (34) of a bank of perfect or near 
perfect reconstruction filters (44, 46, 48). 

The method of claim 1 wherein the step of analyzing 
(36) the signals in each subband comprises the step 
of calculating (64) the average absolute amplitude 
(39) of the signals in each subband and determining 
if any signal peak (41 ) absolute amplitude within the 
subband exceeds the average amplitude within that 
subband by more than a predetermined noise 
threshold amount. 

The method of claim 4 wherein said predetermined 
noise threshold amount is programmable. 

The method of claim 1 wherein said transmitted sig- 
nal is a code division multiple access signal. 

The method of claim 1 wherein the step of analyzing 
(36) the signals in each subband comprises the step 
of calculating (64) the average power level of the 
signals in each subband and determining if any sig- 
nal peak absolute amplitude within the subband ex- 
ceeds the average power level within that subband 
by more than a predetermined noise threshold 
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amount. 

8. The method of claim 1 wherein the step of analyzing 
(36) the signals in each subband comprises the step 
of calculating (64) a DFT or FFT of the signals in 
each subband and determining if any frequency 
component within the subband has a magnitude 
which exceeds the average amplitude level within 
that subband by more than a predetermined noise 
threshold amount. 

9. The method of claim 1 wherein the step of analyzing 
(36) the signals in each subband comprises the step 
of calculating (64) the variance of the amplitude of 
the signals in each subband at each frequency from 
the average absolute amplitude or average power 
level within said subband and determining if the var- 
iance from the average of the amplitude or power 
level at any frequency within the subband has a 
magnitude which exceeds an acceptable level of 
variance. 

10. The method of claim 9 wherein the acceptable level 
of variance is programmable. 

11. The method of claim 1 wherein the step of removing 
(36) ail narrowband interference sources from all 
subbands comprises the step of setting (66, 68, 70) 
the amplitude of all signals in any subband with a 
narrowband interference source present to zero. 

1 2. The method of claim 1 wherein the step of removing 
(36) all narrowband interference sources from all 
subbands comprises the step of adapting the coef- 
ficients of a digital narrowband notch filter so as to 
have a center frequency substantially equal to the 
center frequency of the narrowband interferences 
source. 

1 3. The method o1 claim 1 wherein the step of recon- 
structing (40) a transmitted signal from the individ- 
ual subband signals comprises passing the sub- 
band signals after narrowband interferences sig- 
nals have been removed or attenuated through a 
bank of synthesis filters (78, 80, 82) and up convert- 
ers (72, 74, 76) which form the second half of a plu- 
rality of perfect or near perfect reconstructions fil- 
ters, one such filter dedicated to each subband. 

14. A narrowband interference signal excision circuit, 
characterized in that it comprises: 

means for dividing (34) a transmitted spectrum 
into a plurality of narrow subbands each with 
overlapping frequency limits with sufficient 
overlap to eliminate blind spots or spots in said 
transmitted spectrum of substantially lower 
sensitivity to detection of narrowband interfer- 



ence sources; 

means for analyzing (36) each subband for the 
presence of narrowband interference signals; 
means for eliminating or attenuating (36) any 
5 and all narrowband interference signals found 

in any of the subbands; and 
means for reconstructing (40) a transmitted 
spectrum from the output(s) of the means for 
eliminating or attenuating. 

10 

15. A narrowband interference signal excision circuit, 
characterized in that it comprises: 

a bank of perfect or near perfect reconstruction 

IS filters (44, 46, 48, 78, 80, 82, 1 06, 1 1 8) connect- 

ed in parallel to have a shared input (31 ) for re- 
ceiving a transmitted signal and a shared out- 
put (128) at which a reconstructed transmitted 
signal is output, said bank comprised of a plu- 

20 rality of individual perfect or near perfect recon- 

struction filters each having an analysis filter 
(44, 46, 48) with an input coupled to said shared 
input (31 ) and an output and a synthesis filter 
(78, 80, 82) having an input and having an out- 

25 put coupled to said shared output (128), said 

bank of filters having an overall frequency re- 
sponse which covers the spectrum of said 
transmitted signals and which is comprised of 
a plurality of individual overlapping frequency 

30 responses of the individual filters, each having 

a different center frequency; and 
a detection and cancellation matrix (36) com- 
prised of a plurality of individual detection and 
cancellation circuits (56, 58, 60), each having 

35 an input coupled to an output of a different anal- 

ysis filter (44 : 46, 48) and an output (63, 65, 67) 
coupled to an input of a different synthesis filter 
(78, 80, 82), each detection and cancellation 
circuit comprising a first means (66) for setting 

40 the amplitude of the input signal to the synthe- 

sis filter to which the detection and cancellation 
circuit is connected to zero or attenuating some 
or all of the frequencies in the signal from the 
output of said analysis filter in a first state, and, 

45 in a second state, coupling the output of said 

analysis filter to the input of the synthesis filter 
without attenuation, and a second means (64) 
for determining if a narrowband interference 
signal is present in the output (61) from the 

50 analysis filter to which it is connected, and, if 

not, causing (68) said first means to enter said 
second state, and, if so, causing (68) said first 
means to enter said first state. 

55 16. The narrowband interference signal excision circuit 
ol claim 1 5 wherein said perfect or near perfect re- 
constructions filters (44, 46, 48, 78, 80, 82) are co- 
sine modulated type perfect reconstruction filters 
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which are 256 in number, each with 31 2 coefficients. 

17. The narrowband interference signal excision circuit 
of claim 1 5 wherein the perfect reconstruction fitters 
are implemented using a polyphase filter implemen- 
tation making use of the Noble Identity such that the 
transmitted input signal samples are first processed 
by decimators (100, 102, 1 04) so as to be decimat- 
ed by a factor of M and then filtered into individual 
subbands using a matrix (106) of polyphase analy- 
sis filters, with the output signals of said detection 
and cancellation circuits (112) being filtered by a 
matrix (118) of polyphase filters to implement said 
bank of synthesis fitters with the outputs of the poly- 
phase synthesis filters up-converted in sample rate 
by a plurality of up-converters (124, 126, 12g). 

18. The narrowband interference signal excision circuit 
of claim 17 wherein said decimators (100, 102, 104) 
and up-converters (124, 126, 129) take the form of 
commutating switches (130, 132). 

19. The narrowband interference signal excision circuit 
of claim 15 wherein said first means in each detec- 
tion and cancellation circuit is a switch (66) having 
an output (63) coupled to the input of the synthesis 
filter to which said detection and cancellation circuit 
is connected and having two inputs (61, 70), a first 
of which (61 ) is coupled to the output of the analysis 
filter to which the detection and cancellation circuit 
is coupled and the second of which (70) is coupled 
to ground, and wherein said switch functions to cou- 
ple said switch output to said first input when no nar- 
rowband interference signal is detected (64) and 
functions to couple said switch output to said sec- 
ond input when a narrowband interference signal is 
detected (64). 

20. The narrowband interference signal excision circuit 
of claim 15 wherein said second means (64) in each 
detection and cancellation circuit (56, 58, 60) is a 
processor programmed to calculate the average ab- 
solute amplitude (3g) of the signals output from said 
analysis filter at each frequency in the subband out- 
put by said analysis filter and to compare the am- 
plitude of the signal output by said analysis filter at 
each frequency to said average and determine 
whether, at any frequency, the signal amplitude ex- 
ceeds said average by more than a predetermined 
amount. 

21. The narrowband interference signal excision circuit 
of claim 1 5 wherein said second means (64) in each 
detection and cancellation circuit (56, 58, 60) is a 
processor programmed to calculate the average 
power of the signals output from said analysis fitter 
at each frequency in the subband output by said 
analysis filter and to compare the power of the sig- 



nal output by said analysis filter at each frequency 
to said average power and determine whether, at 
any frequency, the signal power exceeds said aver- 
age power by more than a predetermined amount. 

22. The narrowband interference signal excision circuit 
of claim 1 5 wherein said second means (64) in each 
detection and cancellation circuit (56, 58, 60) is a 
processor programmed to calculate the variance of 
the absolute amplitude or power of each of the sig- 
nals output from said analysis filter at each frequen- 
cy in the subband output by said analysis filter and 
to compare the variance of the signal output by said 
analysis filter at each frequency to the average sig- 
nal power or amplitude, as appropriate, and deter- 
mine whether, at any frequency, the signal variance 
exceeds said average by more than a predeter- 
mined amount. 

23. The narrowband interference signal excision circuit 
of claim 1 5 wherein said second means (64) in each 
detection and cancellation circuit (56, 58, 60) is a 
processor programmed to calculate the average of 
the absolute amplitude or power of each of the sig- 
nals output from said analysis filter at each frequen- 
cy in the subband output by said analysis filter and 
to compare the average of the subband to the av- 
erages of all the other subbands, and determine 
whether the bin average is greater than the other 
bin averages by more than a predetermined 
amount. 
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and generate a colored noise cancellation signal to re- 
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